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1 .Introduction 

There has been recently considerable theoretical and experimental interest in the con- 
cept of free electron lasers (FELs) [1,2]. It has been shown that free electron lasers can 
operate due to different radiation processes: "magnetic bremsstrahlung" in the undulator, 
Smith-Purcel and Cherenkov radiations, radiation in the laser wave. Independing on spon- 
taneous radiation mechanism being principle for a definite FEL scheme, all existing FEL 
devices use for the feedback forming either two parallel mirrors placed at the ends of the 
working area or one-dimentional diffraction grating in which transmitted and diffracted 
(reflected) waves propagate along the electron beam velocity direction (one-dimentional 
distributed feedback (DFB)) (Fig. 1,2). 
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Figure 1: Free electron laser. 
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Figure 2: Smith-Purcel FEL (orotron). 
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According [3], the dispersion equation of the FEL in the collective interaction regime 
is reduced to that of the conventional travelling wave amplifier [4] and the FEL gain at the 
conditions of synchronism is proportional to ^y^, where po is the density of the electron 
beam. 

The volume FEL has been suggested as one of the alternative schemes of FEL which 
provides possibility to design compact sources in various spectral ranges including ultra- 
violet and X-ray [5-10]. 

The main peculiarity of VFEL is the use of one, two or three-dimensional grating as 
a volume resonator providing three-dimensional distributed feedback (Fig.3-5). 
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Figure 3: Parametric quasi- Cherenkov VFEL. 

It is important to emphasize that even one-dimensional diffraction grating may provide 
non-one-dimensional (volume) feedback if the diffracted wave moves in nonback direction 
(the Bragg diffractive angle does not equal |). This results in essential modification of 
the VFEL gain and lasing processes providing, under specific conditions, more effective 
radiation process as compared with conventional FELs using one-dimensional distributed 
feedback. The VFEL gain at the conditions of synchronism is proportional to p])^^^ 
where S is the number of diffracted waves. Volume FEL, if realised, could be made with 
much more compact device structure compared with the conventional FEL and therefore, 
may be interesting for applications in different wavelength regions: from submillimeter to 
X-ray [5-9]. 

It should be emphasized that a fast destruction of the synchronism condition between 
a particle and an emitted electromagnetic wave is characteristic for the VFEL scheme 
with an electron beam passing through a diffraction grating. This leads to the essential 
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Figure 4: Parametric quasi- Cherenkov VFEL. 
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Figure 5: Surface VFEL. 
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increase of the generation threshold parameters. The reduction of the influence of multiple 
scattering becomes possible when electron beam moves either in the split of a grating 
(vacuum VFEL) or over a surface of a grating (surface FEL, SFEL) at a distance ci < A7 
(Fig.3-5)(A is the photon wave length, 7 is the Lorentz factor) [10]. The SFEL has been 
studied in [8]. 

It is easy to understand that vacuum VFEL turns into SFEL when width of the grating 
split grows. 

Radiation mechanisms being the basis of VFEL and SFEL can be various (Cherenkov, 
Smith-Purcel and so on). The spontaneous surface parametric radiation (SSPR) [10] may 
be used for SFEL, for example. We should distinguish the SSPR from Smith-Purcel 
radiation [11,12]. 

The difference between these two types of radiation may be shown by analysing the 
radiation frequency dependence on electron energy. In the case of Smith-Purcel radiation 
the photon frequency is proportional to: 
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and, for photons emitted at small angles to the electron velocity, the radiation frequency 
depends on the electron energy as 7^. Moreover, this wave propagates in vacuum. 
In the case of SSPR, the frequency of photons emitted even at small angles to the electron 
velocity does not practically depend on the electron energy but is determined by the 
Bragg condition. This radiation propagates inside the grating and leaves it only through 
a grating- vacuum boundary. The microscopic nature of both types of radiation is similar: 
they are stipulated by the medium atoms polarisation caused by an electromagnetic fleld 
of a moving charged particle. 
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Figure 6: Smith-Purcel FEL and parametric quasi-Cherenkov VFEL comparison. 



In the Smith-Purcel FEL (orotron, diffraction radiation generator) (see [13-17] ) an 
electron beam passes over a reflecting diffraction grating, two mirrors (or diffraction grat- 
ing) are used for one-dimensional feedback forming (Fig. 2). 
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In the VFEL the non-one-dimensional feedback forming by the diffraction grating is 
used (Fig. 3-5). 

In the present paper the equations describing the VFEL lasing in case of an electron 
beam moving either in a split of a diffraction grating or in a vacuum waveguide containing 
a diffraction grating (vacuum VFEL)(Fig.3,4) have been obtained. The dispersion equa- 
tion allowing to find the vacuum VFEL gain in one-mode generation regime have been 
considered. 

2. Basic formulas describing vacuum VFEL lasing. 

The interaction of an electron beam and an electromagnetic wave propagating along 
a waveguide in the vacuum VFEL is described by the Maxwell and electron movement 
equations 

^ 4 

curl curl E(r,uj) -e{f,uj)E{r,uj) = — —j(r,uj) (1) 

c c 

div e{r,uj)E{f,uj) — 47rp(r, a;) (2) 
— iu!p{f, cu) + divj{f, cu) — (3) 

— * . — * — * 

where E{f,uj) — J e^^E{r,t)dt is the Fourier transformation of the electric field E{f,t); 
s{f,uj) is the dielectric susceptibility of the diffraction grating; j{f,uj) and p{f,uj) are the 
Fourier transformations of the electric current density j{f,t) and electric charge density 
of the beam p{f,t), respectively. 

jif,t)^eY^Vait)S{f-fa{t)) (4) 

a 

p{f,t)^ej:5if-fUt)) (5) 

a 

fa,Va(t) are electron radius- vector and velocity. The subscript a denotes the particle's 
number. 

Movement equations can be written in the form 



dvait) e { ft,,-, , ^ 1 



dt mj I c 



— <j E{fa{t),t) + - Va (t) X H{{f^{t),t) - -^{VaEifait),t) } , (6) 



— * — * 

where E{fa{t),t) and H{fa{t),t) are the electric field and the magnetic field of the elec- 
tromagnetic wave in the point fa(t) at the time moment t, 7 = (1 — ■^)~2. 

Let us consider a sheet electron beam passed over a diffraction grating placed in a 
plane waveguide. At first view this generator is similar to Smith-Purcel FEL (orotron 
or diffraction radiation generator) [3, 13-16]. But, in the volume FEL the radiated wave 
wavelength A is of the same order as the diffraction grating period, the wave undergoes 
Bragg diffraction on the Bragg angle non-equal to | and the diffraction grating provides 
the volume distribution feedback. Let the {y, z) coordinate plane be parallel to the waveg- 
uide (diffraction grating) surface. In the absence of the electron beam the current j = 
and the density p = 0. Equations (1, 2) become periodic in y,z directions (they are not 
periodic in x-direction) . In this case the waveguide dielectric susceptibility can be written 
as: 

e (f, u) = £0 {x) + X (r, u) , (7) 
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where (x) = 1 in vacuum, Eq (x) = Eq in the area of the grating disposition, x ^) is 
the space periodic permittivity discribing the diffraction grating. 
Let the permittivity x ^) be the periodic function of y and z: 

where f/ = ye2 + ze^, is the two-dimensional vector, 62(3) is the unit vector along y{z) axis 
and r = Tye2 + t^Cs is the reciprocal lattice vector of the diffraction grating. 

In the case x = the set of equations (1, 2) desribes passing of the electromagnetic 
waves in the plain waveguide which contains the layer of the matter. The dielectric 
susceptebility of the matter is Eq. The waveguide eigenmodes F„ (x)^and eigenvalues Kn 
are well known [18]. They can be used for simplifying the three-dimensional Maxwell 
equations (1). 

First of all, let us rewrite equation (1) as 

- A ^ (f, cu) - V ( V ((^^0 (x) - 1) E (f, a;)) ) - V ( V (x (r) E (f 



2 22 
-^^ (£0 {x) -l)E (f) - "^x {r) E (f) - "^E (f, u) = (8) 



Equations (8,6) allow us to find the electromagnetic field E{f,u) radiated by electron 
beam. It is well known that gain and the generation threshold can be find in the linear 
approximation. In this case the beam current j is the linear function of E{f, uj): j = jo+Sj, 
where jo is the beam current not perturbated by the radiated field, 6j ~ Eif^oj) is 
the beam current induced by the radiated field. In the linear approximation the set 
of movement equations (6) may be solved by the following way: the electromagnetic 
field E {ra{t),uj) in the right side of equations (6) can be represented as E (r^o + ut,uj), 
where fa (t) — ^ao + ut, fao is the initial coordinate of the electron, u is the electron 
velocity;^^ (t) ~ u in the absence of radiated field. 

As a result, we can obtain from (6) that 

Sg^(uj) = ^ I ^e''^'''-^l—^E(k',uj + k'u] + (9) 

+ (^^-^] (uE(k',u; + k' 
\uj + k'u c^J ^ ^ 

6fa i^) = —SVa (uj) . (10) 
U 

The beam current induced by the radiated field is 

e-'^V"*(5i(r,t)rfVdcu = (11) 




6j{k,uj^ = 

= e ^ g-i'^^t.o I^Sva{uJ - ku) - iu (kSfa{uj - ku)^ | . 

a 

After substitution of expressions (9, 10, 11) in equation (|) we shall obtain the set of 
equations for the field E (r, oj). 
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Let us accomplish the Fourier transformation of the field E {f,uj): 



(12) 



To obtain one-dimensional equation for the field E{x, k\\) let us substitute expansion (|T^) 
in equation (||). 

Let X = (the smooth waveguide) and the electron beam is absent. Then, equation 
(I) allows us to find eigenfuctions ^^(x, k\\) and eigenvalues K^(fc||): 



92 ^ - d 
■— Y'„(x, fc||) - ei — 
ox^ ox 



-ik\\ 



d_ r 
dx 
d_ 

dx 



{eo{x) - 1) Yn,j,{x, fc||)J + i {eQ{x) - 1) A;||K„(x, k\\) 
{£o{x) - 1) K,(x, A^ii) + i {eo{x) - 1) A^||F„(x, A;||) 



- (13) 



^(50 {x)-l)Yn = Kl (A;||)f„. 



If the vacuum-matter boundary is sharp, expression (0) gives the well known equation 
for the waveguide containing the dielectric layer: 
a) in vacuum 

Y„ = kIY„,, (14) 



b) in medium 



Q2 ^ 

-Q^i^n - (^0 - 1) Yn = 



(15) 



Now we can decompose the field E{x,k\^) in terms of the waveguides eigenf unctions 
Yn(x, k\\): 



E {x,k\\) = Y.c^{k\\)\Y^{x,k\\)) . 



The field -^(r', uj) may be represented as 

1 



E(r,uj] 



(2n) 



Cn{k\\) F„(x,A;||))e*'=ii'?rf2fc 



(16) 



(17) 



Let us substitute decomposition (16) into (8) and study the right side of (8) (which is 
determined by the current j{f,uj) = jo{f,uj) + 6j{f,uj)) more attentively. The set of 
equations (8) is a linear system. As a result, we can omit the nonperturbative part of 
current jo and study (8) containing the induced current 6j only. Decomposition (16) 
allows us to obtain the following expression for the right side of (8): 



M 



Airiuj 



Yn{x, k\\ 



-iknfj 



6j (f, u) + -^V (V(5j(f, Lj)) I dxd^ 



Airiuj 1 
c2 (27r)3 



f„(x,fc||)|e*^^ii'? 5j(^,^) - ^k{k6j{k,uj)) 



e'^^'d'^kdx 



;i8) 



A-aiuj 1 
c2 2^ 



j j {Yn{x,k\\)\[5]{k,u) - ^k{M]{k,u)^^ 
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where k = fcy j . It should be mentioned that the electron beam current density j(r, u) 
is not equal to zero only in the vacuum area restricted by the beam transverse size /. This 
fact results in the appearance of integral / in (18) which means the integration over the 



area where j{f,uj) ^ (see Fig.3,4). As a result, we have 



M 



Aniuj 1 
c2 2^ 



dkro (Yn{k^,k\\) I 6j{k,u;) - ^k (k5j{k,uj)) \ , 



(19) 



where 



Yn{kx,k\\) — / (Yn{kx,k 



The expression for the current density 6j{k,u;) contains the sum F — ^O'^ao 

a 

US average this sum over distribution of the particles in the beam: 



(20) 



u u 

where ^{kx — k'^) — J e~'^^^''~^'=''>^i^{x)dx, y / (p{x)dx = 1, Hs the characteristic transver- 

a a 

sal size of the beam, function (p{x) describes the distribution of the particles along x- 
direction, no is the electron density of the beam. 
Using (20) we can write 5j{k,ijj) as 



Sj{k,u!) = udip{k,u!) = u 



dx'e Lp{x') X 
(uE{x',k\\,uj)) . 



(21) 



As a result, using (21) we can represent (19) as 



M 



Amu 1 
c2 2^ 



dkx (Yn{kx, fcy) < u - —^K^^ \ ^^{k, i^) 



(22) 



Attiu} 1 



c2 27r 



^2 ^ 



u -k > 6ip{k^ uj) 



dkx (^En (j^xj k 

In the cold beam case (when the condition kxUx\ <S L to be fulfilled) we can write 
Amuj ie^UQ 1 



M 



(ix (F„(a;, /cii, a;) \u k\\—i— e\ \ x (23) 

\ \ Ci; " ox UJ 



(d d c^oj"^ 1 \ 



where operator ^ acts on the functions disposed on its left. 
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Substituting the decomposition E{x', k^, uj) =J2 Cn'(^ll) Yn{x, kn)) into (23) and using 

n ' 

the ortogonahty of the eigenf unctions we obtain from (8, 23) 



k^-(^- K 



T,n' 



c 



(24) 



ie noAr, 



n ' (to — k\\itj rwycu 



Cn'ihl 



where 



dx 



f c^- d 

Yn{x,kii,uj) <u ^11 -^(t^ ei) 

uj ox uj 



X 



d ^^^^ + 
dx dx M^7^ 



(m|f^(x, k\\,uj) 



27rc2 



dx 



X 



Yn{x) 



u k\\ 

UJ 



X (25) 



where Xe//( ll'^ll ''") effective permittivity. It contains two terms: the first pro- 

portional to x(^) ^) and the second proportional to V (Vx(^, ^) 



h + r) 



Xr[.X) 



r„/(x,/cii +f)) + 



(26) 



dx{Yn{x,k\\) —kikxAx) F„/(x, ^n + f) 



~* ft ~* 
where k = + ik\i. 

Let the condition k\\ S> ^ be fuUfiled {d is the diametrical size of the waveguide). In 
this case the second term of (26) is less than the first one. 

We can simplify the system (2) solving by assuming a practically important case of a 
single mode n existance in the waveguide. It is possible when condition a;^/c^x""''(fc||, fcy + 
r) (k^ — f^n') ^ ^ 1 is fulfilled. In this case all terms with n' 7^ n in the sum on n' in (24) 
can be omitted. 

As a result, we obtain the set of equations which is similar to that describing the 
multiwave dynamical diffraction of electromagnetic waves in a diffraction grating. 

In the two wave diffraction case the Bragg condition accomplishes for two waves with 
wavevectors k\\ and : fc'i ~ + , 



and the set of equations (|^) can be 



written as 



UJ 



k\\ + 



ujlA, 



7C^ (uj — k\\u 



UJ 



Cnik\\) - —Xeffih^h+^Mk\\ + T) = 0, 



(27) 



UJ 



-x:Uh + ^^h)(^n{k\\) + 



UJ 



2n{k\\ + f) = 0, 



where Eq = 1 — c k^/uj\ ujl is the Lengmuer frequency of the electron beam {uj\ 
Ane'^nQ/m). 
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The set of equations (27) is similar to that for the electromagnetic field amplitudes 
describing lasing in volume FEL for the case of beam moving in volume diffraction grating 
[6] . The main discrepancy appeares in the dependence of the equations (27) on k\\ . The 
similar set of equations for volume FEL depends on k. 

As a result, we can conclude that all the main results obtained for VFELs holds true 
for the vacuum VFEL. First of all, the non-trivial solution of system (27) exists only when 
the determinant of the system is equal to zero. This allows us to obtain the dispersion 
equation for fcy and a;: 

{co - k\\uy [(fcjc^ - uho) {{h +t)c'- uho) - uj\rxZ] = (28) 



According to [5-8] the study of the dispersion equation let us find the condition of the 
appearance of the convection and absolute instability of the beam and, as a result, obtain 
the gain and the generation threshold. 

In the plain waveguide case (x™ = 0) from (28) we have 



2 

kiiuY (^/cjc^ - = -—Ann (29) 



UJ 
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(u - klluY (fc||C - UJ^) ~ Ann (30) 

Equations (^), (|30|) coincide with the equations describing the wave spectrum for 
Cherenkov instability of the beam in medium. From (|30D we can obtain that k\\ has 
imaginary part Imk\\ and 



2 



" 2c V 2cjeo7 / 

when the Cherenkov condition 1 — lo j c^fsQCOs'Q = is fulfilled. As we see, lmk\\ is 

1 /3 

proportional to ng (no is the density of the electron beam). It means that the gain is 
proportional to Uq"^ as well. This dependence is typical for all types of one-dimensional 
FEL in the collective regime [3]. 

Let the waveguide contains a diffraction grating (x"*^ ^ 0). The wave spectrum is 
described by equation (^). When the coefficient Ann = 0, equation (|2^) splites into two 
equations: 

(kY - + f)' - - uj^Txr, = D {^,u;) = (32) 



(uj-kuy = (33) 



Equation (^) describes the electromagnetic wave spectrum for the waiveguide containing 
a diffraction grating. Equation (|33|) describes the wave spectrum of the electron beam 
charge density. Let us study the solutions of (|28|) near the point where the left side of 
I p^) is equal to zero. The solution of (31, 32) in the vicinity of the exact Bragg condition 
k\\ + T ~ k\\ can be written in the form 

fc,o = (1 + (5) , cuo = fc°M(l + <5), (34) 
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where S <^ 1 and can be found from exact Bragg conditions 

= _?fe±l!, (35) 
z axis is parallel to the beam velocity u. From (28,34) we can obtain for S 



Now we can study equations (g^)- Let us rewrite (|3) as 

{uj — kzuf D {kz,uj) = A {kz,uj) , (37) 
where ^ 

A {k,, UJ) = - Y^n" ((^11 + r)' _ ^2^^^ ^3g) 

The solution of ( |57| ) can be represented as 

c<j = cjq + c^', = kzo + k'^ (39) 

where \uj'\ <^ ujq and <^ |A;^o|- 
Let us write D {kz, uj) in the form 



It can be shown that (dD/du)^^ ^ ^ can not be equal to zero. That is why, we can omit the 
term proportional to {d"^ D / duJ^) ^ On the other hand, the derivative {dD/dkz) ^ 
can be equal to zero. In this case equation (36) can be written as 

[uj - k^u) [k^ -Fuj) = , (41) 

afc? 



where 

fdD\ fd^DV^ 

F = -2 



du I \ dkl I 



For uj' equation (40) takes form 



u \ dkl 
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As a result, equation (|4^) gives 



Imk., 



^0 \Q\ \Xt\ r 



1/4 



1/4 



2r2 



(43) 



Let us remind that the Langmuer frequency ul ~ "n-o • So, accordingly to (^31), we 
have been obtained a very important result: in vacuum VFEL both Imkz and the gain are 
proportional to rig''^ for the definite orientation of the diffraction grating in the waveguide 
in contrast with the conventional one-dimensional FEL for which the gain is proportional 
to n^^. From (32) and (43) we have: 



Imkz 
Imkii 



> 1, 



(44) 



because 



Ujl\A„n\ 



■C 1, ^ ^ 1 and Ann -C • As a result, in our case of the volume 



feedback the gain is larger then that for one-dimensional feedback. For example, the de- 
pendence of the threshold current density for the volume and one-dimensional geometries 
on the length is represented in Fig. 7. 
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Figure 7: The dependence of the threshold current density for the volume (1) and 

o 

one-dimensional (2) geometries on the length (A = 6283 A)- 



The gain becomes higher when the distributed feedback is formed by the multi-wave 
dynamical diffraction: 



where 5* is the number of diffracted waves. 



uj /wl, |A„„| V+i 



(45) 



12 



CONCLUSION 



The vacuum VFEL amplification and the generation process develop more intensively 
then in ordinary FEL using one-dimensional distributed feedback. Such the VFEL, if 

realised, could be made with much more compact device structure compared with the 
FEL and therefore, may be interesting for applications in different wavelength ranges: 
from submillimeter to X-ray. Such the VFEL can be realised on the basis of the existing 
accelerators. 

Author is gratefuU to K. Batrakov for the help in the threshold current densities 
computer simulation carrying out. 
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